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The brain is much less hard-wired as traditionally thought. 
Permanently, new synapses are formed, existing synapses 
are deleted or connectivity rewires by re-routing axonal 
branches (structural plasticity). However, all current large- 
scale neuronal network models are hard-wired with plasti- 
city merely arising from changes in the strength of existing 
synapses, therefore missing an important aspect of the 
plasticity of brain networks. This project is to develop the 
first large-scale neuronal network model with structural 
plasticity in the neuronal network simulator NEST [1] and 
to make it scalable for HPC. 

Formation and deletion of synapses in the model for 
structural plasticity (MSP) [2] depends on the number of 



synaptic contact possibilities that each neuron has, i.e. the 
number of axonal boutons and dendritic spines. Therefore, 
we developed a framework that allows the addition of 
synaptic elements (i.e. axonal boutons or dendritic spines) 
for every neuron model already implemented in NEST. 
The user can then define its own synaptic elements and 
their corresponding growth dynamic depending on the 
electrical activity (see Figure 1). Synapses are formed by 
merging corresponding synaptic elements or are deleted 
when synaptic elements are lost. The update in connectiv- 
ity depends on the availability of the synaptic elements in 
the entire networks. To make this model scalable for HPC, 
we developed a probabilistic approach that reduce both 



import nest 

# Create one neuron 

neuron = nGst.create( "iaf_neiiron" ) 

# Dictionary to describe the growth dynajnic 

# of the synaptic elements 
Dynamic - { 

# growth curve type 

' growth_curve ' : "gaussian", 

# initial nurijer of synaptic elements 
' z_niinus ' : 0.0, 

# scaling of the growth curve 
■ growth_rate ' : 0.001, 

# Parameters of the growth curve 
■eta': 0.1, 

■eps': 0.7, 

) 

# Add two synaptic elements using the previous 

# tiynamlc to the neuron 
nest. Set Status (neuron, 

■ synaptic_elements ' , 
{'Axon': gc, 'Dendrite': go}) 
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Figure 1 Example of the definition of two synaptic elements for an integrate and fire neuron using PyNEST. Here, the growth curve of the 
synaptic element is a Gaussian defined by three parameters: u (the growth rate), r|z and e (intersections with the x-axis). The growth dynamic 
depends on the electrical activity of the neuron modeled by the intracellular concentration of calcium. 
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communication between compute nodes and their mem- 
ory usage. 

This implementation of the MSP in NEST allows neu- 
roscientists to address important scientific questions on 
how large-scale networks rewire their connectivity in 
response to distortions in electrical activity balances 
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